Determining the crystal-field ground state in rare earth Heavy Fermion materials 

using soft-x-ray absorption spectroscopy 
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We infer that soft-x-ray absorption spectroscopy is a versatile method for the determination of the 
crystal-field ground state symmetry of rare earth Heavy Fermion systems, complementing neutron 
scattering. Using realistic and universal parameters, we provide a theoretical mapping between the 
polarization dependence of Ce M4,5 spectra and the charge distribution of the Ce 4/ states. The 
experimental resolution can be orders of magnitude larger than the 4/ crystal field splitting itself. To 
demonstrate the experimental feasibility of the method, we investigated CePd2Si2, thereby settling 
an existing disagreement about its crystal-field ground state. 

PACS numbers: 71.27.-l-a, 75.10.Dg, 75.30.Mb, 78.70.Dm 



Heavy Fermion materials are strongly correlated rare 
earth or actinide materials where the atomic-like / elec- 
trons interact with conduction electrons giving rise to 
extraordinary low energy properties. There has been a 
revival of interest in many of these well studied com- 
pounds with an intense search for new ones since the 
discovery that many of these materials are at the bor- 
der between magnetic order and superconductivity when 
driven through a quantum critical point by either ap- 
plying a magnetic field or pressure, and/or exhibit non- 
Fermi liquid behavior [J, i, i, i, [sj, [ej, [tI, i, Q . There are 
indications that the unconventional superconductivity in 
these compounds is triggered by antifcrromagnetic cor- 
relations leading to the interesting aspect of competing 
magnetic and superconducting order parameters which 
may involve a momentum dependent hybridization of / 
and conduction electrons 13, 11, l3]- Therefore, any mi- 
croscopic understanding of how the wealth of properties 
in Heavy Fermion materials evolves out of the strongly 
interacting / ground state requires the knowledge of the 
spatial distribution of the / state involved. It has been 
pointed out already early by Zwicknagl, for example, that 
the Fermi surface of several s yste ms depends strongly on 
the crystal-field symmetries [13|. Other effects like the 
quadrupolar orderingreflect the charge distribution of 
the / ground state 14, [l^ [H, 17 1, but so far rehable 
input from experiment is scarce. 

The standard experimental technique to determine the 
crystal-field energy level scheme is inelastic neutron scat- 
tering on polycrystalline samples. The additional infor- 
mation of quasi- and inelastic intensities should be suf- 
ficient to extract information about the / spatial distri- 
bution as well, but the analysis of magnetic intensities is 
often hampered by broadened lines, phonons in the same 



energy window as the magnetic excitations or strong ab- 
sorption of one of the sample's constituents (e.g. Rh, In, 
B, or Sm). Therefore, the wave functions of the crystal- 
field states are often determined from a combined analy- 
sis of neutron and single crystal static susceptibility data 
illli,!!!. However, aspointed out by Witte et al. ^ 
and Janouasova et al. [22], even this combined method 
often fails due to the powder and thermal averaging in 
either technique. Moreover, in the presence of magnetic 
order anisotropic molecular field parameters have to be 
introduced in order to fit the static susceptibility, ending 
up with too many free parameters for a unique descrip- 
tion of the crystal-field potential, leaving ambiguities and 
unresolved debates in the present literature. One then 
has to resort to time consuming inelastic polarized neu- 
tron scattering experiments with their demand for large 
single crystals [2l|, |22|, or, equally time consuming, to 
elastic scattering experiments with polarized neutrons to 
determine the / magnetic form factor [2^. The latter 
has the advantage that it is also applicable for systems in 



which crystal-field excitations are not defined [2j|. Nev- 
ertheless, neutron based methods have their limitations 
in case the compound under study contains more than 
one type of magnetic ion. 

We will show in this report at the example of CePd2Si2 
that polarization dependent soft-x-ray absorption spec- 
troscopy at the rare earth M4_5 edges (soft-XAS) is a 
powerful tool to give undisguised direct information con- 
cerning the 4/ charge distribution in the ground state. 
This technique is complementary to inelastic and elas- 
tic neutron scattering, but has the advantage of having 
orders of magnitude better signal-to-noise/background- 
ratio and of requiring only small amounts of sample ma- 
terial. It is capable of measuring materials which are 
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FIG. 1: (color online) Calculated M4,5 soft-x-ray absorption 
spectra for pure \Jz) states of Ce^"*" for the incoming light 
polarized || and _L c (z || c in tetragonal symmetry - see below). 
On the right the spatial distribution of the 4/ electrons for 
the respective | J^) states. 



strongly absorbing for neutrons, and it is also element 
specific, so that the presence of more than one type of 
magnetic ion is no limitation. 

Over the last 20 years soft-XAS has developed into a 
powerful technique for studying the element-specific lo- 
cal electronic structure of 3d transition metal and 4/ rare 
earth compounds. The respective 2p-to-3(i (L2,3)and 3d- 
to-4/ {M4_5) absorption processes can be well understood 
in terms of atomic-like transitions into multiplet-split fi- 
nal states 



26l . l27l . |28|. Unique is that the dipole 



selection rules are very effective in determining which of 
those final states can be reached and with what intensity 
for a given initial state symmetry. This makes the tech- 
nique an extremely sensitive local probe, ideal to study 
the valence, spin and orbital state of the ions of interest. 

The majority of the soft-XAS studies on rare earth ma- 
terials were so far focused on issues concerning their mag- 
netism by measuring the magnetic circular dichroic effect 
[2^. Surprisingly, only little attention had been paid on 
orbital symmetry or crystal-field issues in rare earths: 
there are a few soft-XAS reports early on addressing sur- 
face and interface crystal-fields 13, 31, l^]- As far as bulk 
materials are concerned, only the giant crystal-field sys - 
tem CeRh3B2 has been investigated 33, 34, H, H, 37 1, 
thereby giving the impression that it would be difficult to 
study experimentally systems with smaller crystal-fields. 
Actually, CeRh3B2 is a quite special case since the large 
size of the crystal-field splitting is comparable to that of 
the spin-orbit interaction, resulting in an intermixing of 
the J=7/2 and J=5/2 orbital states [38|. Here we fol- 
low up on the theoretical work of Jo et al. [nil and 
show that for the small crystal-fields in rare earth sys- 
tems a simple mapping can be found between the ground 
state orbital symmetry and the polarization dependence 
of the spectra, allowing an accurate and quick quantita- 



tive analysis. 

The spin orbit interaction in the / shell of rare earth 
compounds is much larger than the crystal-field splitting 
due to its small spatial distribution. J in the LS coupling 
is therefore a good quantum number, and the eigenstates 
in the presence of crystal-fields can be expressed in terms 
of linear combinations of the corresponding | J^) states. 
In the case of Ce'^"'" with in tetragonal site symmetry, 
for example, the 6-fold degenerate J=5/2 state splits into 
three doublets: |1) = a|±5/2) - 5|t3/2), |2) = |±l/2) 
and |3) = 6|±5/2) + a|T3/2) with |r^), \Tq), and \T'^) 
symmetry, respectively. The crucial question is which 
one of these states forms the ground state and what are 
the a and h coefficients (a^-f6^=l), thus defining the ac- 
tual charge distribution. In Stevens formalism [s^ the 
tetragonal crystal-field potential of an is given by the 
three crystal-field parameters B2, and Bj, and can 
be expressed alternatively through the transition ener- 
gies AE12 and Ai?i3, and the coefficient a. 

To study the sensitivity of the Ce M4^5 spectra to the 
4/ charge distribution of the ground state, we performed 
atomic calculations which include the full multiplet the- 
26 
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21, \M, \2M using the XTLS 8.3 programm 



It accounts for the intra-atomic 4/-4/ and 3d-4/ 
Coulomb and exchange interactions, the 3d and 4/ spin- 
orbit couplings, and the local crystal-field parameters. 
All atomic parameters were given by the Hartree-Fock 
values, with a reduction of about 35% for the 4/-4/ pa- 
rameters and 25% for the 3d-4/ to reproduce best the 
experimental isotropic spectra of Ce materials with low 
Kondo temperatures, as to account for configuration in- 
teraction effects not included in the Hartree-Fock scheme 

We start with calculating the M4.5 spectra for the pure 
\Jz) states of the Ce'^+ with 4/^ and J=5/2 for two po- 
larizations of the electric field vector E, namely || c 
and E ± c, where the tetragonal c-axis is aligned along 
the z direction. Fig. 1 shows the results of these calcula- 
tions. The spectra are dominated by the Ce 3d core- hole 
spin-orbit coupling which splits the spectrum roughly in 
two parts, namely the M5 {hiy « 877 - 887 eV) and M4 
{hi/ « 895 — 905 eV) white lines regions. Crucial is that 
there is a strong polarization dependence which is char- 
acteristic for each of the states. The existence of such 
a polarization dependence is intuitively clear when look- 
ing at the spatial distributions of the / electrons of the 
respective \Jz) states as shown also in Fig.l: the distribu- 
tions differ remarkably between the different \Jz) states. 

The results in Fig.l are a nice demonstration of the 
so-called initial state symmetry effect, an important con- 
cept for the determination of the ground state orbital 
symmetry in real materials: the splitting between the 
various crystal-field states can be small, e.g. smaller than 
the broadening of the spectra due to lifetime, phonon or 
experimental resolution, and yet the polarization depen- 



3 




890 895 900 905 910 

Photon energy in eV 

FIG. 2: (color online) Top: Calculated M4,5 soft-x-ray ab- 
sorption spectra for the mixed states a|±5/2) - ti|=p3/2) for 
mixing factors from to 1 for light polarized || and _L c. 



dence can be very large. Important is that the tempera- 
ture is sufficiently low so that primarily only the lowest 
state is populated and contributes to the signal. In fact, 
we can make use of the fact that the crystal-field split- 
tings in rare earth materials are typically much smaller 
than the «400 meV spectral broadening. When calcu- 
lating the polarization dependence of the states |1) or 
3), we find they can be well approximated by an inco- 
herent sum of the |3/2) and |5/2) spectra, weighted with 
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FIG. 3; (color online) Bottom: experimental Mi,^ soft-x-ray 
absorption spectra of CePd2Si2 at T = 50 K for light polarized 
II and 1. c. Top: simulated M4,5 spectra with a? =0.29. 



the proper a? and 6^ parameters. The results are shown 
in Fig. 2, where we let a? run from to 1, with a^=0 
corresponding to a pure |3/2) and a^ = l to a pure |5/2) 
state. Fig. 2 in effect can serve as a map to find the 
coefficient directly from an experiment. It should be 
generally applicable for tetragonal Ce materials without 
the need to know the crystal-field splittings themselves as 
long as the total crystal-field splitting is small in compar- 
ison to the spin orbit splitting of 230 meV. Yet, having 
knowledge of the crystal-field energies e.g. from neutron 
scattering, this XAS technique will also enable us to give 
a clear and unique set of crystal-field parameters, inde- 
pendent of molecular fields. Important is, of course, that 
this atomic approach is only accurate for materials in 
which the Kondo temperature is sufficiently smaller than 
the splitting between the lowest and first excited crystal- 
field state, otherwise the analysis needs to be extended 
by using the Anderson impurity model. 

We have chosen to measure the linear polarized soft- 
XAS spectra of CePd2Si2 in order to proof the feasi- 
bility of the technique. Several groups measured the 
crystal-fields scheme of CePd2Si2 with neutron scattering 
[S 41 1 and agree in the energy splitting but end 

up with considerably different ground state wave func- 
tions. The authors of and [l^] both took the energy 



splitting from neutron data and fitted the single crystal 
static susceptibility in order to obtain the ground state 
wave function. Both groups obtained decent fits although 
different ground state coefficients were assumed. This 
is due to the fact that anisotropic molecular field pa- 
rameters were introduced in order to take into account 
the antiferromagnetic order at 10 K, ending up with too 
many free parameters for a unique fit. The CePd2Si2 
single crystal used was grown by the Czochralski method 
as described in [11]. It orders at Tn « 10 K. The XAS 
measurements were performed at the Dragon beamline of 
the NSRRC in Taiwan. The spectra were recorded using 
the total electron yield method in an ultra-high- vacuum 
chamber with a base pressure of 2xl0~^" mbar. Clean 
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FIG. 4: (color online) Temperature dependence of the exper- 
imental M4,5 soft-x-ray absorption spectra of CePd2Si2. The 
inset gives the crystal-field transition energies as in [40^ 



sample areas were obtained by cleaving the crystals in- 
situ. The photon energy resolution at the Ce M4.5 edges 
{hv « 875 — 910 eV) was set at 0.4 eV, and the degree 
of linear polarization was k, 98%. The CePd2Si2 sin- 
gle crystal was mounted with the c-axis perpendicular to 
the Poynting vector of the light. By rotating the sample 
around this Poynting vector, the polarization of the elec- 
tric field vector can be varied continuously from E \\ c to 
E ^- c. This measurement geometry allows for an optical 
path of the incoming beam which is independent of the 
polarization, guaranteeing a reliable comparison of the 
spectral line shapes as a function of polarization. 

The bottom curves of Fig. 3 show the measured soft- 
XAS spectra of CePd2Si2 at T = 50 K for the light polar- 
ized parallel and perpendicular to the c-axis. The spec- 
tral features are well resolved and their sharpness can be 
taken as a confirmation for the atomic character of the 
Ce states. Important is that there is a polarization de- 
pendence and that this can be clearly seen due to the 
excellent signal-to-noise/background ratio. Comparing 
these experimental spectra to the simulations in Figs. 1 
and 2, one can immediately deduce that the ground state 
must be state |1) with a? close to 0.3. Indeed, a more de- 
tailed simulation analysis yields the best fit for a^=0.29. 
We know from neutron data [i^, |4l| that the first excited 
crystal-field level is above 18 meV, i.e. well above 200 K, 
so that at T = 50 K only the ground state contributes to 
the spectra. When simulating the data the Boltzmann 
factor was nevertheless taken into account. The value 
found for a? corresponds to a=0.54 and 6=0.84, so that 
we essentially confirm the ground state of CePd2Si2 as 
proposed by 0: |1) = 0.55 |±5/2) - 0.84 |±3/2). The 
corresponding spatial distribution of the 4/ electron in 



this ground state is shown in the inset of Fig. 3. 

To confirm the crystal-field nature of the polarization 
effect, we have also carried out measurements at elevated 
temperatures. The polarization dependence in the spec- 
tra diminishes as the temperature rises as is shown in Fig. 
4 for 50, 100, 150, 200 and 295 K. This is in accordance 
with the increased population of higher crystal-field lev- 
els, resulting eventually in an isotropic spectrum at suf- 
ficiently high temperatures. These elevated temperature 
spectra can be well described using the ground state sym- 
metry and coefficient from the 50 K spectrum and the 
crystal-field energies from neutron data [i^, depicted in 
the inset of Fig. 4. One only has to take into account 
the Boltzmann factor for each respective temperature. 

In conclusion we have shown that the Ce M4_5 spec- 
tra have a polarization dependence which is characteris- 
tic for the charge distribution of the Ce 4/ states. We 
have provided a set of reference spectra (Figs. 1 and 2) 
from which one can find the ground state parameters di- 
rectly from the experiment, without the need to know 
the crystal-field energies themselves as long as they are 
small. Using this method we have investigated the local 
electronic structure CePd2Si2 and settled the debate con- 
cerning the ground state wave function of this compound, 
demonstrating that soft-x-ray absorption spectroscopy is 
a powerful method complementing neutron scattering. 

We would like to thank Lucie Hamdan for her skillful 
assistance in preparing the experiments, and John My- 
dosh and Daniel Khomskii for valuable discussions. 
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